Conical flow due to partonic jets in central Au+Au collisions 
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In jet quenching, a hard QCD parton, before fragmenting into a jet of hadrons, deposits a fraction 
of its energy in the medium, leading to suppressed production of high-pT hadrons. The process can 
generate shock waves. We study the distortion of Mach shock waves due to jet quenching in central 
Au+Au collisions and its effect on particle production. Finite fluid velocity and inhomogeneity of the 
medium can distort the Mach shock front significantly such that the inside shock front disappear and 
the outside shock front is opened up. We also show that the STAR data on azimuthal distribution of 
background subtracted secondaries, associated with high pr trigger, are reasonably well explained 
by the excess pions produced due to partonic energy loss. 

PACS numbers: PACS numbers: 25.75.-q, 13.85.Hd, 13.87.-a 



The three most important results that came out from 
the heavy ion programme at RHIC are (i) dramatic 
suppression of inclusive hadrons production at large 
transverse momentum (high p? suppression) p], [3, [ij], 
(ii)disappearance of away side two hadron correlation 
peak Q and (iii) large elliptic flow [B| . High px suppres- 
sion confirmed the theoretical prediction of jet quench- 
ing 0] . Long before the RHIC Au+Au collisions, it was 
predicted, in a pQCD calculation, that in a deconfincd 
medium, high-speed partons will suffer energy loss, lead- 
ing to suppressed production of hadrons. The observed 
high pt suppression in Au+Au collisions are in agree- 
ment with the prediction. Large elliptic flow observed in 
non-central Au+Au collisions confirms fluid like behavior 
of the produced matter. The elliptic flow data are well 
explained in an ideal hydrodynamic model, with initial 
energy density of deconfined matter £j ~ 30 GeV/ fm 3 , 
thermalized at an initial time Ti=0.6 fm Q. Transport 
based models, on the other hand, require unrealistically 
large cross-section to reproduce the observed elliptic flow 
Q . All these observations are being treated as evidences 
for the creation of a very dense, color opaque medium of 
deconfined quarks and gluons [6[. 

As the jet quenching models suggest, if the partons 
lose energy in the medium, a natural question arises, 
what happens to the lost energy? In [9|, it has been 
suggested that a fraction of the lost energy goes to col- 
lective excitation, called "conical flow". The quenching 
jet moves at speed of light (cj et w 1), much larger than 
the speed of sound of the medium (c 2 s < 1). The quench- 
ing jet can produce a shock wave with Mach cone angle, 
6m = cos~ 1 c s /cj e t- The resulting conical flow will have 
characteristic peaks at <fi = tt — 9m and <fi = tt + 9m- Indi- 
cations of such peaks are seen in azimuthal distribution 
of secondaries associated with high px trigger in central 
Au+Au collisions 13, U|- As Mach cone is sensitive to 
the speed of sound of the medium, it raises the possibility 
of measuring the speed of sound of the deconfined matter 
of Quark-Gluon-Plasma. 

Recently, we have (numerically) solved hydrodynami- 



cal equations with an additional source, representing the 
quenching jet [13]. We concluded that for realistic jet 
energy loss, the detailed structure seen in the azimuthal 
correlation of secondaries, can not be explained by Mach 
cone formation. However, that simulation was limited 
to b=0 collision. More importantly, only one trajectory 
(along the diagonal) was considered. Recently, distortion 
of Mach regions, due to radial and longitudinal expansion 
was studied [131 ]. Finite fluid velocity induces significant 
deformation of Mach shocks. For the perturbation (jets) 
moving along the diagonal, the Mach cone opening angle 
(6m = sin~ 1 c s /cj e t = 7r/2 — 6m) is reduced, but it re- 
mains symmetric, 9m+ =0m-i where, 9m+ and 9m- are 
the upper and lower cone opening angles, with respect 
to the jet. Mach cone become asymmetric when the jet 
moves along a chord (9m+ ^ 9m- ) ■ The relation between 
9m+ and 6m- reversed as the jet moves from upper half 
plane to lower half plane. It is important see, whether 
exact numerical simulation corroborate such predictions. 
Such changes in Mach angle, will also affect the particle 
production as opposed to the only diagonal jets studied 
in [13]. 

A schematic representation of the jet moving through 
the medium is shown in Fig[TJ We assume that just be- 
fore hydrodynamics become applicable, a pair of high-py 
partons is produced at P. One of jet moves outward and 
escapes, forming the trigger jet. The other enters into 
the fireball. For jet trajectory parallel to the x-axis, the 
production point can be characterised by the angle <fip ro d 
varying between [— tt/2, +tt/2]. Strong jet quenching and 
survival of the trigger jet, forbids production in the inte- 
rior of the fireball. Jet pairs can be produced only on a 
thin shell on the surface of the fireball. For an impact pa- 
rameter b collision, we assume that the pair is produced 
on the surface of the ellipsoidal fireball with minor and 
major axis, A = R - 6/2 and B = R^Jl - b 2 /4R 2 with 
R=6.4 fm (for Au+Au collisions). The fireball is expand- 
ing and cooling. The ingoing parton travels at the speed 
of light and loses energy in the fireball which thermalizes 
and acts as a source of energy and momentum for the 
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fireball medium. 

The model we use is described in [l2| ■ Briefly, we solve 
the energy-momentum conservation equation, 



where the source is modeled as, 



J v {x) = J{x) (1,-1,0,0), 
S 3 (r 
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Massless partons have light-like 4-momentum, so the cur- 
rent J v describing the 4-momentum lost and deposited 
in the medium by the fast parton is taken to be light- 
like, too. rjet(t) is the trajectory of the jet moving with 
speed \dx- ]e t/dt \ —c. jr-{x) is the energy loss rate of the 
parton as it moves through the liquid. It depends on the 
fluid's local rest frame particle density. Taking guidance 
from the phenomenological analysis of parton energy loss 
observed in Au+Au collisions at RHIC [14| we take 
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where s(x) is the local entropy density without the jet. 
The measured suppression of high-p^ particle production 
in Au+Au collisions at RHIC was shown to be consistent 
with a parton energy loss of | = 14 GeV/fm at a ref- 
erence entropy density of sq = 140 fm -3 We will use 
this realistic jet energy loss through out the present cal- 
culation. 

For the hydrodynamic evolution we use a modified ver- 
sion of the publicly available hydrodynamic code AZHY- 
DRO @, El- The code is formulated in (r,x,y,r}) co- 
ordinates, where r—\/t 2 ~z 2 is the longitudinal proper 

time, ?7=^ln |^| is space-time rapidity, and rj_ = (x, y) 

defines the plane transverse to the beam direction z. 
AZHYDRO employs longitudinal boost invariance along 
z but this is violated by the source term We there- 
fore modify the latter by replacing the ^-function in ([3]) 
by 
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with a = 0.70 fm. 

Intuitively, this replaces the "needle" (jet) pushing 
through the medium at one point by a "knife" cutting 
the medium along its entire length along the beam direc- 
tion. Thus assumption of boost-invariance will over esti- 
mate the effect of jet quenching. While a complete study 



of this would require a full (3+l)-dimensional hydrody- 
namic calculation, the present boost-invariant simulation 
should give a robust upper limit for the effect produced 
by the quenching jet. 

The modified hydrodynamic equations in (r, x, y, rf) co- 
ordinates read [H, HBj 

d T T TT +d X (V X T TT )+dy(iyf TT ) = -P+J, (6) 
d T T TX + d X (V X T TX ) + d y (V y f TX ) = -d X P~J, (7) 

drf^ + d x (v x f T *) + d y (Vyf Ty ) = -d y p, (8) 

where f^ = rT^, i % =T Ti /T TT , p = rp, and J = tJ. 

Y 




FIG. 1: Schematic representation of a jet moving through 
the medium. The high pr pair is assumed to produce on the 
surface of the fireball characterized by the angle 4> pr od- One 
of the jet escapes forming the trigger jet, the other enters into 
the fireball. 

We solve these equations for Au+Au collisions at im- 
pact parameter b — 3 fm. Roughly it corresponds to 
0-5% centrality Au+Au collisions. We use the stan- 
dard initialization described in 0] and provided in the 
downloaded AZHYDRO input file [HI], corresponding 
to a peak initial energy density of eo = 30 GeV/fm 3 at 
To = 0.6fm/c. We use the equation of state EOS-Q de- 
scribed in 0, HH incorporating a first order phase tran- 
sition and hadronic chemical freeze-out at a critical tem- 
perature T c = 164 MeV. The hadronic sector of EOS-Q is 
soft with a squared speed of sound c 2 s w 0.15. 

In Figl^l we have shown the results of our simulations. 
We have considered all positions of the jet trajectories, 
— it/2 < tfiprod < 7r/2. Some representative cases, jets in 
the upper half plane {§ pro d = 45°), in lower half plane 
{<j)prod — —45°) and along the diagonal, are shown in 
Figl21 In the left panels we have shown the constant en- 
ergy density contours after r=8.6 fm of evolution. With 
realistic energy loss, the perturbation to energy density 
is small and conical flow developed by the Mach shock 
wave is barely visible. To see the Mach cone formation, 
we subtract the energy density of the fluid evolved under 
a similar condition, but with out the quenching jet. Con- 
tours of excess energy density, due to the quenching jet, 
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FIG. 2: Left panels (a), (b) and(c) shows the constant energy 
density contours at time t—8.6 fm, The jet is produced at 
4> V rod = —45°, 0° and 45°. The right panels, (d), (e) and 
(f ) show the excess energy density due to quenching jet only. 
The conical structure is clearly visible. Approximate Mach 
cones are depicted by the solid lines. To show the distortion, 
for the off diagonal jets, we have shown the symmetric Mach 
cone by the dotted lines. 



are shown in the right panels, d, e and f. Mach cone like 
structure is clearly visible. We also note that the distor- 
tion of the Mach cone due to transverse velocity roughly 
corresponds to the theoretical study [13| ■ In the figure, 
we have drawn the (approximate) Mach cone. For the 
jet moving along the diagonal, the Mach region is sym- 
metric with respect to the jet axis, the upper and lower 
cone opening angles are same, 0m+ = @M— ■ But for the 
jets moving along a chord, it is asymmetrical. Due to 
transverse expansion Mach region is pushed out. The in- 
ner surface is pushed more than the outer surface. In the 
upper half plane, 6m+ < @m-, the relation is reversed in 
the lower plane, 6m+ > 9m -■ While the simulations cor- 
roborate the predicted distortions of Mach shock regions, 
we also note that the Mach region is more complex than 
predicted in Mach surface is not smooth. Inhomo- 
geneity of the medium plays significant role is distorting 
the Mach region. 

Using the standard Cooper-Frey prescription, we now 
calculate px integrated pion spectra 4? , at the freeze-out 
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FIG. 3: Azimuthal distribution of pions due to quenching jet 
alone. In the upper panel, the quenching jet is confined to 
upper half plane, 4> pr od~0° to 75°, in steps of 15° (from top 
to bottom). In the lower panel the jet is confined to lower 
half plane, 4> pro d= —75° to 0°in steps of 15°. The arrows 
indicate the Mach peaks (n ± 9m) for resonance hadron gas 
(c 2 s a 0.15). 



surface at freeze-out temperature of 120 MeV. In Figj3^i, 
for various position of the jet (confined to upper half 
plane) angular distribution of the pions is shown. To 
show clearly the effect of deposited energy, we have sub- 
tracted the distribution obtained in a evolution without 
a quenching jet. When the jet is at the periphery of 
the upper half plane, 4> pro d = 75°, the angular distri- 
bution is hardly changed in the lower half plane. Par- 
ticle production is enhanced at ^ w 0.85rad. We also 
observe depleted production at <j> = 4> pro d- The distribu- 
tion donot show any second peak. As the jet trajectory 
comes closer to the centre, qualitatively, the distribution 
remains unchanged. The peak position remains nearly 
unaltered but it gets widened significantly. The deple- 
tion of production around <fi = <f> pro d continues. We also 
find that as the trajectories comes closer and closer to the 
diagonal, particle distributions in the lower half plane is 
also increased. However, the second peak eludes us. For 
the trajectories closer to the diagonal, the single peak 
observed in off diagonal jet trajectories, gives way to a 
broad maxima. A mirror image is found for the jets in 
the lower half plane, Fig[3jD. The jet enhances the par- 
ticle production around <j) « 4rad. As before, the peak 
gives way to a broad maxima, as the jet trajectory comes 
closer and closer to the diagonal. Here also, we find de- 
pleted production at <p = 4>prod- Apparently, as the jet 
moves in, it drags fluid along with it, resulting a reduced 
particle production around <f> = 4> pro d- 

In the Fig[3l the arrows indicate the Mack peak po- 
sitions, cf> — 7r ± 8m, for a static resonance hadron gas 
with squared speed of sound cj, ~ 0.15. If we asso- 
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ciate the excess production with shock wave formation 
due to quenching jet [l6j |. interesting conclusions can be 
drawn. For a jet confined to upper half plane, the shock 
wave is distorted to such an extent that Om- ~ and 
&M+ ~ 0.72rad, while for the jets confined to lower half 
plane, 6m+ ~ and 9m- ~ 0.72rad. The angles should 
be compared with the cone opening angle 6m± ~ 0.4rari 
for resonance hadron gas. Due to finite fluid velocity 
and inhomgeneity of the medium, the inside shock front 
disappear, while the outside shock front opened up sig- 
nificantly. 
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FIG. 4: Black circles are STAR data on background sub- 
tracted secondaries in azimuthal angle A<f>, associated with 
hard trigger (pr > 4 GeV), in 0-5% centrality Au+Au colli- 
sions. The solid line is the normalised azimuthal distribution 
of excess pions, due to quenching jet only, in b = 3 fm Au+Au 
collisions. 

As mentioned in the beginning, STAR collaboration, 
in 5% central Au+Au collisions, has seen structures in 
the azimuthal distribution of background subtracted sec- 
ondaries, associated with high px {pt > 4 GeV) trigger 
[l(| ■ Are such structures associated with conical flow as 
suggested in Q? In Figf4]we have compared the prelimi- 
nary STAR data (Io| with (normalised) azimuthal distri- 
bution of pions due to quenching jet only. The spectrum 
is obtained after averaging over all the possible jet tra- 
jectories. The away side azimuthal distribution of back 
ground subtracted STAR data closely match with the an- 
gular distribution of pions due to quenching jet only. We 
also reproduce the hint of depression at <ft = ir. 

To summarise, we have studied the Mach shock wave 
formation in central Au+Au collisions due to a quenching 
jet. Hydrodynamical equations, with a source (represent- 
ing the quenching jet) are numerically solved with initial 
conditions appropriate for Au+Au collisions at b=3 fm 
and with jet energy loss, consistent with observed high 
Pt suppression. Explicit simulations indicate that the 
Mach shock front generated by the quenching jet, depend 
significantly on the jet trajectory. Finite fluid velocity 
(and also inhomogeneity of the medium) cause signifi- 



cant distortions of the shock region. For a jet moving off- 
diagonally, Mach cone become asymmetric with respect 
to jet axis. The inside shock front is pushed in, while the 
outside shock front is pushed out. Azimuthal distribu- 
tion of pions, due to quenching jet only, suggest that for 
a off-diagonal jet trajectory, the distortions can be large, 
such that the inside Mach shock front disappear. We 
have also shown that the (normalised) azimuthal distri- 
bution of pions due to quenching jets only, averaged over 
all the trajectories, reasonably well explains the STAR 
data [l(| on the background subtracted secondaries, as- 
sociated with high pt trigger, along with the hint of de- 
pression at <f> = 7T- It strongly suggest that the structure 
is closely linked with "conical" flow due to partonic en- 
ergy loss. 
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medium is changed. 



